Double gap and solitonic excitations in the spin-Peierls chain CuGe03 



M. Ain 1 *, J.E. Lorenzo 2 , L. P. Regnault 3 , G. Dhalenne 4 , 
A. Revcolevschi 4 , B. Hennion 1 and Th. Jolicoeur 5 . 

1 Laboratoire Leon Brillouin (CEA-CNRS) CE-Saclay. 91191 Gif/Yvette Cedex FRANCE. 
2 European Synchroton Radiation Facility, 38043 Grenoble FRANCE. 
3 Laboratoire de Magnetisme et de Diffraction Neutronique, CENG. 38054 Grenoble Cedex 9 
FRANCE. 

4 Laboratoire de Chimie des Solides (Bat 414) - UA CNRS 44@ Universite de Paris-Sud. 91405 
Orsay Cedex FRANCE. 

5 Service de Physique Theorique, CE-Saclay. 91191 Gif/Yvette Cedex FRANCE. 



We have studied magnetic excitations in the dimerized spin-Peierls phase of 
CuGe03, by high resolution inelastic neutron scattering. We measured the well- 
defined spin triplet dispersive mode which is gapped throughout the whole Brillouin 
zone. We also observed that this mode is separated by an unexpected second gap of 
~ 2 meV from a continuum of magnetic excitations extending to higher energy. The 
first gap (or 'triplet gap') and its associated dispersive mode is due to the break- 
ing of a singlet dimer into a delocalized triplet. We propose that the second gap 
(or 'solitonic gap') and the continuum correspond to dissociation of that triplet into 
two unbound spin- 1/2 solitons that are separated by a dimerized region of arbitrary 
length. 



PACS numbers: 75.10.Jm, 75.40. Gb, 64.70.Kb 



Typeset Using REVTEX 



1 



The recent observation by Hase |]] of a characteristic magnetic susceptibility in CuGeOs, 
dropping abruptly to zero below T 5 p=14.3 K, clearly suggested that it was a new one- 
dimensional spin-Peierls compound. This was confirmed by X-Rays photographs []2]|| 
that revealed, below Tsp, superlattice peaks indexing according to the propagation vec- 
tor ksp= (0.5,0,0.5). These two experimental evidences indicate : 1- That a gap has opened 
over a non magnetic singlet ground state as demonstrated by neutron studies |4|,[| and 2- 
That the crystal was undergoing a magnetoelastic distortion where copper ions dimerize 
with their left or right nearest neighbor along the chains. As a result, the initially uniform 
exchange coupling becomes staggered. 

Single crystals of CuGe03 are grown in an image furnace by the traveling floating zone 
method. They belong to the orthorhombic space group Pbmm. Magnetic chains of Cu ++ , 
S=l/2 ions are parallel to the c axis. The spin-Peierls gap is observable at ku,F=(0,0,0.5) 
or equivalent points, but there are no magnetic Bragg peaks. Dispersion curves of magnetic 
excitations along the three principal directions a,b and c are of simple sinusoidal shape [|]||; 
estimates for intrachain nearest neighbors exchange (NNE) along c gave J\ ~ 120 K, as 
derived from fits to classical magnon theory in ref or to Bonner-Blote relation in ref || . 
Fits to classical magnon theory |4],||, indicate that NNE between chains are only one and 
two orders of magnitude smaller along the b and a directions respectively. 

In this letter we present experimental evidence that there are in fact two gaps in this 
system and not only one as predicted by the classical approach . Although this observa- 
tion could fit in the framework of Cross-Fisher theory ||, we shall use a solitonic approach 
to elaborate an explanation which accounts for two gaps. A dimerized system has an obvi- 
ous excitation which consists of breaking a dimer bond into a triplet at a cost of a certain 
magnetoelastic energy. The triplet will be delocalized along the chain generating eigenstates 
of definite momentum. However there is another possible excitation in this system because 
the triplet can absorb a second amount of energy corresponding to a second gap and thus 
dissociate into two S=l/2 traveling solitons that generate the continuum. This has some 
analogies to the well known two-spinon continuum of the uniform (undimerized) Heisenberg 
S=l/2 AF chain (HAFC) that has been investigated extensively in KCuF 3 (see figure 6 of 
ref 0). Above T$p the analogy is even more complete as we shall see on figure |3]. On the 
other hand it has been suggested that competing next-nearest-neighbor exchange (NNNE) 
was the driving mechanism in the dimerization process of CuGe03 instead of magnetoelas- 
tic coupling as generally admitted. In support of this idea, were satisfactory fits ||T0Hl2| of 
the magnetic susceptibility of CuGe03 above Tsp- Since this susceptibility was not well 



reproduced by the Bonner-Fisher curve |13 , which is appropriate to the isolated S=l/2 AF 
NNE chain, incorporation of NNNE gave better results. Yet in the final part of section II 
of our discussion we express some reservation with respect to this interpretation based on 
competing NNNE. 

Inelastic neutron scattering measurements have been performed on two triple-axis spec- 
trometers : 1- 4F1 (Orphee reactor, LLB Saclay) has an incident beam, fixed in direction, 
extracted by a pair of graphite monochromators (the second one being vertically focussing). 
It was operated at constant k/=1.55 A (5.01 meV) with a horizontally focussing graphite 
analyser and a berylium filter to cut out higher-order components of the diffracted beam. 2- 
IN14 (HFR, ILL Grenoble) has one vertically focussing graphite monochromator and a hori- 
zontally focussing graphite analyzer, the rest of the set up was similar on both spectrometers. 



2 



The resolution on both apparatus was very close to 0.2 meV (FWHM) as deduced from the 
incoherent peak at zero energy transfer. In both experiments the same single crystal (nearly 
1 cm 3 ) was oriented with the b and c crystallographic axes in the scattering plane. Two 
series of inelastic scans were recorded for neutron energy transfers ranging from -0.3 meV 
to 11.5 meV. All scans are corrected for A/2 contamination in the incident beam 

Figure ([!]) shows three energy scans at T=2.6 K. They correspond to excitations near 
the zone boundary along the c* direction parallel to the chains. Five elements are visible 
on the scan at Q=(0,l,0.5) : 1- The zero-energy incoherent peak showing the spectrometer 
resolution. 2- A first gap, called hereafter the 'triplet gap' with a value of A = 2 meV at 
Q = (0, 1,0.5). 3- A well-defined magnon-like mode first observed by Nishi et ai. Q]. This 
mode is in fact a spin-triplet mode as shown by measurements in a magnetic field ||. Its 
asymmetric shape is due to convolution of instrumental resolution and steep curvature of the 
dispersion curve in the vicinity of Q=(0,l,0.5). 4- The intensity between the middle peak 
(or triplet mode) and the plateau, falls to the background level. This is clearly a new gap in 
energy that we call hereafter the 'solitonic gap'. At Q=(0,l,0.5), this 'solitonic gap' is close 
to 2 meV. Defining the background in this experiment, is an issue that will be addressed 
when presenting figure @. 5- Finally, brought out by the 'solitonic gap' we find some 
intensity (low but clearly present) which constitute the expected continuum that extends at 
least up to the maximum energy transfer of our study, i.e. 11.5 meV. Scans for Q=(0,l,0.48) 
and Q=(0,l,0.46) display the same structure as that for Q=(0,l,0.5) with the exception of 



the incoherent peak which is not shown. We recall that neutron |5]U and Ramann [15 
scattering have already given clear evidence for the existence of such a continuum. 

Figure (Q) shows a series of six energy scans regularly spaced along b* between 
Q=(0,l,0.5) and Q=(0,2,0.5) at T=1.7 K. No incoherent peak here, only the peak of the 
dispersive mode followed again by the 'solitonic gap' and the continuum are visible in this 
series of scans. Note that owing to coupling between chains, as already mentionned, there is 
dispersion along and the positions of the peaks are not constant in energy as they would 
be for a pure one-dimensional system. 

Figure (H) provides a more detailed picture of the energy scans at Q=(0,l,0.5). It shows 
that when the temperature is raised to T=29 K the peak of the dispersive mode drops 
and widens, filling in the double gap region and merging with the continuum indicating 
that the two gaps have collapsed and that we have recovered the continuum of the uniform 
HAFC. When we reach T=150 K, the magnon-like mode and the continuum have totally 
disappeared, we are then in the truly paramagnetic region; note that the intensity falls to 
the level of what was measured at T=1.7 K in either gap, this level is considered as the 
background of our experiment. In the insert we subtracted the scan at 150 K from the scan 
at 1.7 K, both were corrected for Bose factor after subtraction of a background of 15 counts 
on each; what remains is the dispersive mode, the 'solitonic gap' and the continuum. A 
phonon is distinguishable near 11 meV in the 150 K data. 

The fact that the double gap has been overlooked in previous experiments []4|-§] is due to 
poorer resolution of the instruments used before. In these former experiments the high energy 
tail of the uncompletely resolved dispersive mode precluded observation of the 'solitonic gap' 
by causing a smooth crossover to the high energy continuum. In the present experiment 
high resolution was obtained through the use of a small kj=1.55 A. 

It has been impossible to detect acoustic phonon branches around k$p, and moreover, 
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preliminary polarized neutron measurements at Q=(0,l,0.5) and Q=(0.5,5,0.5) indicate that 
all of the intensity in the peak of the dispersive mode is magnetic, as is the major part if not 
all, of the intensity in the continuum. The absence of an inelastic nuclear contribution is 
consistent with the fact that in the dimerized phase of CuGeC>3, intensities of nuclear super- 
lattice peaks at {ksp} are very weak and displacements of atoms extremely tiny. Nonetheless 
superlattice nuclear peaks are sensitive to a magnetic field as proved by the commensurate- 
incommensurate transition that occurs at 12.5 T |T7[. All this suggests that there could 



exist a magnetoelastic coupling through spin-charge hybridization jlfj . 

There are three models that are related to our observations in CuGeC>3 : I - The AF 
chain with NNE J x and NNNE J 2 . II - The AF chain with NNE only but with imposed 
staggering of the exchange. Ill - The AF chain coupled to a phonon field u(x) |§. Although 
all three models predict a gap, only model III predicts a double gap which is consistent 
with our experimental results on the excitation spectrum of CuGeO"3. An analysis of these 
models has been given by Haldane [ITH . 



I - In the case of NNE J\ and NNNE J 2 , the Hamiltonian of the chain has full translational 
invariance. If J 2 is smaller than a critical coupling (J2/ J\ < 0.2412(1)), the ground state is a 
spin liquid but if J 2 is larger, the ground state is dimerized and twice degenerate. Translation 
invariance by one lattice spacing is spontaneously broken. We will refer to this situation as 
"spontaneous dimers". The effective long- wavelength, low-energy theory is described by a 



sine-Gordon model 19 



H = J dx^{U 2 + (V0) 2 ) + acos(/30). (1) 

In this equation, II is the momentum conjugate to <fi which is related to the z-component of 
the spin at position x by : S z {x) = — V0(x)/v / 27r + C(— ) x cos(v / 27T0(x)). When J 2 7^ , the 
value of the sine-Gordon coupling is (3 = 2y/2n. If J 2 is larger than the critical value, one 
is in the massive (or gapped) phase of the theory of Eq. (JJ). The (3 = 2v / 2~7r sine-Gordon 
theory has no bound states [^(| an d the elementary excitations are kinks that correspond 
to a ±27T variation of the argument of the cosine term in Eq. ([I]) over a localized region 
of space. These solitons therefore have spin S=l/2. The physical picture is simple : an 
excitation means that a singlet in the dimerized ground state is broken into a triplet that 
immediately disintegrates into two free solitons. As a consequence the magnetic excitations 
form a continuum above some threshold and there is no well-defined mode below. This is 
not what we observe in our experiments. 

II - We turn now to the externally dimerized chain: there is an additional modulated 
exchange SJ2 n (~) n S n • Sn+i- There we have explicit doubling of the unit cell. This also 
leads to a sine-Gordon model Eq. (0) but now with a coupling (3 = y/2n. The kinks that 
still correspond to a ±2tt variation of the argument of the cosine in Eq. ([!]) have now spin 
S z = ±1. The sine-Gordon theory with (3 = \[2tx has two breather bound states [^0j : 



one which is degenerate with the kink states. This state completes the S=l triplet which 
is expected due to the full rotational invariance of the theory. The other bound state is a 
singlet and thus plays no role in the magnetic excitation spectrum. Here the physical picture 
is quite different from that of model I. The singlet bonds are pinned to the lattice by the 
dimerizing potential. The elementary excitation corresponds to breaking a singlet bond in 
a triplet and then this triplet will move along the chain. This triplet state is not a domain 
wall and cannot disintegrate as has been seen in numerical simulations [ 101,111 J2~I |. There 



4 



are continua above this well-defined mode that are due to excitations of several triplets. 
Haas and Dagotto have recently performed a study of the dynamical properties of an 
externally dimerized chain including a NNNE J 2 . They have shown that there is a continuum 
starting immediately above the spin triplet mode contrary to our finding of a 'solitonic gap' 
in CuGeOa. 

To the extent that NNNE should be visible on the shape of the dispersion curve it 
becomes informative to calculate the dynamics of a colinear S=l/2 Heisenberg AF model 
with J i < and J 2 , and lattice spacing c. It yields the following dispersion relation : 



hu(q) = 2| singc|y (Ji(Ji - 4J 2 ) + 4J 2 sin 2 gc (2) 

which is only valid for J2/J1 < 0.25 (Villain's criterion). We see that the curve u(q) is 
narrowed by the last term with the J 2 factor, which makes it depart markedly from a 
sinusoidal shape in contrast to our observations on CuGeC>3. The analogy with CuGeOs is 
enforced by noting that the underlying AF system in CuGeO"3 would have AF Bragg peaks 
at kAp=(0,0,0.5) as verified in CuGeo.993Sio.oo703 |22|] where AF and dimerization coexist. 



(If J 2 > and J 2 /Ji > 0.25 the calculation ought to be conducted differently as appropriate 
for helimagnetic incommensurate AF). 

Ill - Finally we turn to the more realistic model having an elastic displacement field 



|8|j23|]. It involves an additional coupling of the sine-Gordon field with an elastic field, of the 
form u(x) cos(v^27T0(a;)). The Hamiltonian Eq.(|l]) again has full translational invariance. 
This invariance is spontaneously broken and there are thus domain walls: the displacement 
u(x) in a domain wall goes from +uq on one side of the chain to — uq on the other side, or 
vice- versa. The spin soliton involves only a variation of 7r of the argument of the cosine and 
thus has spin S=l/2, as in the case of the spontaneously dimerized chain. This elementary 
soliton can be visualized as an isolated S=l/2 copper spin that separates two regions of the 
chain that are dimerized. These solitons have been studied in the past [ 23[| . No detailed 



information is available on the sine-Gordon theory coupled to an additional scalar field, 
however, since the coupling is well in the massive regime, it is clear that we expect bound 
states of these solitons. The most likely candidate is a triplet of the same nature as in the 
externally dimerized chain. If enough energy is available to overcome the binding energy, 
this triplet state can then disintegrate into two solitons. We expect then a triplet mode that 
is well-defined below the solitonic continuum. This is consistent with what we observe. 

To summarize, we have shown by inelastic neutron measurement that there is a mid- 
gap dispersive mode and confirmed the existence of a continuum of excitations. We have 
proposed that this continuum is made of unbound S=l/2 domain walls. 

It is a pleasure to thank S. Aubry, M. Azzouz, A. R. Bishop, G. J. Mclntyre, M. Poirier, 
H. Schulz and A. Tsvelik for interesting discussions and M. Geoghegan for reading the 
manuscript. 
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FIGURES 



FIG. 1. (IN14-ILL) Three energy scans for Q=(0,1,Q C ) with Q c ={0.5, 0.48, 0.46}, at T=2.6 
K. They are vertically shifted apart for clarity. The horizontal graduation is common to all scans; 
the left vertical axis is for Q=(0,l,0.5) scan only. Each horizontal arrow indicate the zero intensity 
level for the scan above it. In the insert a general view of the first scan displaying the five elements 
described in the text. Now labelling a dimer in its singlet state by • — •, the triplet state by ff and 
a spin 1/2 on a copper site by f, we can represent the peak of the magnon-like mode as a traveling 
triplet • — • • — • f|~ • — • • — •, then after the 'solitonic gap' the continuum would correspond 
to delocalized spins 1/2 such as • — • j • — • • — • f • — •. 

FIG. 2. (4F1-LLB) Six energy scans for Q=(0,Q 6 ,0.5) with Q 6 ={1, 1.2, 1.4, 1.6, 1.8, 2} at 
T=1.7 K. Same convention for axes as on Fig (|l]). Maximum peak position is indicated by a 
vertical arrow, the intensity reached is written below. The sharp peak of the magnon-like mode, 
the 'solitonic gap' and the continuum are clearly visible in all scans. 

FIG. 3. (4F1-LLB) Three energy scans at Q=(0,l,0.5). 1- At 1.7 K (circles) we have succes- 
sively, the incoherent peak, the 'triplet gap', the magnon-like mode that reaches 1983 counts at 
2 meV, the 'solitonic gap', and the continuum. 2- At 29 K (diamonds) the continuum, similar to 
that of the S=l/2 HAFC. 3- At 150 K (squares) the purely paramagnetic region. In the insert, 
subtraction of the scan at 1.7 K from the one at 150 K showing the magnon-like mode and the 
continuum. 
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